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We provide precise constraints on the size of any black holes forming in the early Universe for a 
variety of formation scenarios. In particular, we prove that the size of the apparent horizon of a 
primordial black hole formed by causal processes in a flat Friedmann universe is considerably smaller 
than the cosmological apparent horizon size for an equation of state p — kp (1/3 < fc < 1). This 
also applies for a stiff equation of state (fc = 1) or for a massless scalar field. The apparent horizon 
of a primordial black hole formed through hydrodynamical processes is also considerably smaller 
. . . than the cosmological apparent horizon for < fc < 1. We derive an expression for the maximum 

size which an overdense region can have without being a separate closed universe rather than part 

■ of our own. Newtonian argument shows that a black hole smaller than the cosmological horizon can 

' never accrete much. 
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I. INTRODUCTION 

In the early stage of the universe, black holes may have formed from cosmological density perturbations J,]. Such 
"primordial" black holes (PBHs) may produce 7-ray background radiation and cosmic rays through quantum evap- 
oration or contribute to the density of the universe, and this leads to observational constraints on their number 
density . Studying the formation and evolution of PBHs therefore leads to important constraints on models of the 
] early universe. However, to calculate these constraints, one needs to know how large PBHs were at their formation 
T-H ■ epoch and how much they subsequently increased their mass through accretion. 

] A simple physical argument suggests that the size of a PBH will initially be comparable with the cosmological 
particle horizon scale for a hard equation of state 0. This is because an overdense region at maximum expansion 
must be larger than the Jeans scale in order to collapse against the pressure but not so large that it forms a separate 
closed universe. Since both these scales are comparable to the particle horizon size, with the first being only slightly 
i-S^ • below the second, the size at maximum expansion must be finely tuned if a PBH is to form. (For a soft equation of 
state, the Jeans scale is much smaller than the cosmological horizon scale, but the collapse is likely to lead to pancake 
Q ■ formation and fragmentation unless the region is very spherically symmetric 0|, so one again needs fine-tuning.) Since 
I the size of a region at maximum expansion is determined by its initial overdensity (compared to the flat Friedmann 
c/J . background), only density fluctuations with a narrow range of amplitudes can produce PBHs. 

The limits on the size of a PBH have only been estimated rather roughly in previous papers. However, the fraction 
of the Universe going into PBHs is very sensitive to these limits 0, especially if the fluctuations on a given scale 
have a Gaussian distribution, so one would like to calculate them more precisely. In order to do this, we first need 
^ . to determine the various cosmological scales that arise in the early Universe: in particular, the cosmological particle 
' horizon, the cosmological apparent horizon, the cosmological sonic horizon and the Jeans length. Depending on the 
formation scenario, these may all be relevant to PBH production. They all have the same order of magnitude but 
they also depend on the equation of state parameter k (where p — kp) and one purpose of this paper is to calculate 
these dependencies for an exact Friedmann model very precisely. 

The region undergoing collapse will not itself be part of a flat Friedmann model, so in order to discuss the size 
of a PBH, one needs to make some assumption about the form of the overdensity. One approach is to assume that 
the universe is exactly Friedmann beyond some matching radius. In this case, the black hole apparent horizon is 
necessarily within the matching radius but its precise position (relative to some cosmological scale) depends on the 
density profile in the overdense region. The simplest models have constant overdensity, in which case the region can 
be regarded as part of a closed Friedmann model. However, more complicated possibilities can be envisaged, so it is 
important to obtain more general results. 

In this paper we will derive upper limits on the black hole size which are independent of any assumption about 
the density profile. In particular, we prove that the black hole apparent horizon size is considerably smaller than 
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the cosmological apparent horizon scale for an equation of state parameter in the range 1/3 < A; < 1. An important 
tool in proving these results is the use of double-null coordinates. This is partly because, with the highly non-linear 
fluctuations required for PBH formation, it is not always clear how to construct spatial hypersurfaces of constant 
time. It is also because double-null coordinates are convenient for numerical calculations. Although we do not present 
such calculations here, we do present them in an accompanying paper ^ for the case in which the universe contains 
a scalar field. 

If the fluctuations which generate the PBHs derive from inflation or quantum gravity effects, they may effectively 
extend to infinity (or at least well beyond the usual Friedmann particle horizon) rather than being confined to a finite 
region. In this case, the above results are inapplicable and the only firm upper limit on the size of the PBH comes 
from the separate universe condition. Although this constraint has been discussed in general terms it cannot be 
formulated precisely without specifying the density profile within the overdense region. However, we will show that 
it can be calculated exactly in terms of the equation of state parameter k if the overdense region is assumed to be 
homogeneous. 

The amount of accretion by a PBH is also sensitive to its initial size. A simple Newtonian calculation suggests 
that accretion is unimportant if the PBH is much smaller than the cosmological horizon at formation but that one 
could have self-similar evolution, with the black hole growing at the same rate as the universe, if the size is initially 
comparable with the cosmological horizon scale. The standard formation criterion implies that this is indeed the case, 
which suggests that there could be a large amount of accretion for a hard equation of state. However, self-similar PBH 
solutions have been studied extensively 8, 9, 10, and it has been found that no such solution contains a black 

hole inside an exact flat Friedmann solution for a perfect fluid with < k < 1. This applies even for a stiff fluid with 
A; = 1, although this has been the subject of some controversy These results have been confirmed by numerical 

calculations Indeed simulations suggest that, the harder the equation of state, the more difficult it 

will be for a PBH to form and the smaller the mass accretion rate will be. Also one can show that a relativistic effect 
suppresses the accretion when the PBH has a size very close to that of the cosmological apparent horizon Q . 

II. DOUBLE-NULL FORMULATION OF EINSTEIN'S FIELD EQUATIONS 

If we adopt units in which G = c = 1 and the abstract index notation of |17| , the Einstein equations 

Rab-^gabR=8TTTab (2.1) 

can be written as 

Rab = Snfab, (2.2) 

where 

Tab = Tab — -^gabT, (2.3) 

T EE T"^. (2.4) 

We consider a spherically symmetric system, for which the line element can be written in the form 

ds^ = -a'^{u, v)dudv + r'^{u, v){d9^ + sin^ ed(p'^), (2.5) 

where u and v are advanced and retarded times, respectively. Then, from Eq. H2.2|l . we obtain the following partial 
differential equations: 

^ = 2^^-4^T„„, (2.6) 
r a r 



— -4^T,„ (2.7) 
r a 

= 2^^-4^r™, (2.8) 
r a r 

2 ^ 



r 



-347rT,,. (2.9) 



There are no other independent components of the stress-energy tensor: spherical symmetry requires Tgg — s\v? 9T^^ 
and Tg^ — 0. The initial data for the above equations are most naturally given on two initial null surfaces, u = uq 
and V = vq. On each null surface the geometry is described by the two functions a and r. 
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The Hawking mass m, which is a well-behaved quasi-local mass in spherically symmetric spacetimes, is defined as 



r 



m^-[l + -^]. (2.10) 



Equations H2.6|) - (|2.9|) then imply (see also fl^ l 

87rr2 



a2 



We adopt the double null coordinate system as a natural choice of coordinates because we will be dealing with 
highly nonlinear perturbations in the cosmological background, so it is not always guaranteed that there exists a 
"cosmologically preferred" comoving time slicing. 



III. UPPER LIMIT ON BLACK HOLE APPARENT HORIZON SIZE 

In the following, we focus on an ingoing null surface v = V[j, as shown schematically in Fig. For a given matter 
stress-energy tensor, a can be chosen arbitrarily on this surface and then r is determined by Eq. (|2.t)|) . This implies 
that the form of the function a ori v — vq can be understood as a gauge freedom. 

At the black hole apparent horizon the expansion along the outgoing null congruence vanishes, i.e. = 0, while 
at the cosmological apparent horizon the expansion along the ingoing null congruence vanishes, i.e. — 0. The 
cosmological apparent horizon is distinct from the cosmological particle horizon and defined in terms of the quasi- local 
properties of spacetime. The radius of the cosmological apparent horizon is the Hubble radius in a flat Friedmann 
universe, as shown later. From the definition ()2.1U|I of the Hawking mass, it is clear that r = 2m at these apparent 
horizons. 

Lemma 1 The condition 2m /r < 1 holds if and only if r^ui",v < 0. The condition 2m /r = 1 holds if and only if 
T',uT.v = 0. The condition 2m/ r > 1 holds if and only if r^u'r.v > 0. 

Proof Trivial from Eq. ifTTUI) . □ 

Theorem 1 Suppose there exists uq such that r > 0, r y > and r „ < at {uq,vq). Then if T^u ^ 0, there exists a 
centre Uc (> uq) such that r{uc,VQ) = andr{u,Vo) is a strictly decreasing function of u for uq < u < Uc ■ 

Proof. Equation (|2.t)|l can be rewritten as 

The condition r > at (uo,vo) implies that r is non-zero there, so the left-hand side is negative on the null surface 
V — vq. Since ?■ > and r,u < at {u,vo) = (uo,wo), we have {r,u/a^){u,VQ) < {r^u/<i^)iuo,vo) < for u > uq and 
r > 0. Since a'^{u,vo) is subject to the gauge choice, we can always assume that a'^{u,vo) is bounded from below by 
a positive value. Thus is bounded from above by a negative value. Therefore, there exists Uc (> uq) such that 
r = at (uc, Vq) and such that r{u, vq) is a strictly decreasing function of u for uq < u < Uc. □ 

Lemma 2 Suppose there exists uq such that r > 0, > and r.„ < at (uq,vq). //T„„ > and there exist roots 
A such that 2m/ r ^ \ at (A, uo) for uo < A < Uc, then < and > for uq < u < ubh, where u = ubh is the 
smallest root and corresponds to the black hole apparent horizon. If there is no such root, then > and r.„ < 
for Uq < u < Uc. 

Proof. Trivial from Lemma ns-nd Theorem ^ □ 

Theorem 2 Suppose there exists uq such that r > 0, > and r.„ < at (uo,wo)- IfT^u > and T^v > 0, then 
m{u,VQ) is a decreasing function of u for uq < u < ubh if o, black hole apparent horizon exists at u = ubh > Uq and 
for uq < u < Uc otherwise, where r{uc, Wq) = 0. m{u, vq) is constant if and only if T^u — T„.„ — 0. 

Proof. Since r^u < and > in the relevant interval, equation H2.11|) implies m^u < 0. □ 
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Theorem 3 Suppose that T„„ > and T„„ > and that there exists uq such that r > 0, r.„ > and r.„ < at 
(uo,Vo). If a black hole apparent horizon exists at ubk > uq, then its area radius is less than or equal to 2m(ito, Wq). 
Equality holds if and only if T^u — T^v — 0. 

Proof. From Lemmanand Theorein|21 we have ^(mbhj vq) — 2m(uBH, vq) < 2m{uo, vq)- □ 



Here we comment on the energy conditions assumed above. For a perfect fluid, we have 

T"'' ^ pu^'u'' +p{g'''' + u''u''), (3.f) 

where u", p and p are the four velocity, energy density and pressure, respectively. The condition Tuu > is equivalent 
to requiring p > —p. The condition Tuv > is equivalent to requiring p > p. These two conditions together are 
equivalent to the dominant energy condition p > \p\. If we assume an equation of state of the form p = kp, the 
dominant energy condition requires — 1 < fc < 1 and p > 0. For a massless scalar field we have 

jnah ^ ^,a^,6 _ l^ab^.c^^^^ (3.2) 

and the conditions r„„ > and T^v = automatically hold. Recall that a massless scalar field can be identified with 
a stiff fluid (with the equation of state p = p) if it is vorticity-free and has a timelike gradient |19| . However, we do 
not need this equivalence in the following discussion. 

IV. APPLICATION OF UPPER LIMIT TO A PBH 

The above considerations are completely general, but we now focus on the situation in which one has a black hole 
in a cosmological background. To make the discussion precise, we adopt a "compensated" PBH model, in which the 
black hole apparent horizon is contained within a perturbed region, which is surrounded by an exact flat Friedmann 
solution beyond some (generally finite) radius. Various approaches have been used in this context. One is to model 
the black hole as part of a closed Friedmann model, surrounded by a Schwarzschild void, so that the total mass 
at the edge of the void corresponds to that within a flat Friedmann model at the same radius Another is to 

match the collapsing region to the flat Friedmann background by a shock-wave or a sound-wave, since this allows 
a discontinuity in the density or density gradient . If the PBH forms from density perturbations generated in a 
preceding inflationary phase, then the black hole could have a size larger than the particle horizon, but we do not 
consider such a situation in this section. 

We therefore have various possible situations, depending on the size of the matching radius. These possibilities 
correspond to the following definitions, in which the upper limit on the matching radius is progressively reduced. 

Definition 1 // the matching between the perturbed region containing the PBH and the flat Friedmann region is made 
at some finite radius, then we describe the PBH as "locally produced". 

If the perturbation results from a causal process, it is natural to assume that the perturbed region is inside the particle 
horizon of the surrounding fiat Friedmann region. It should be noted that "causal", in this sense, is distinct from 
"physical" , since there exist physical mechanisms (such as inflation) for producing non-causal perturbations. 

Definition 2 // the matching between the perturbed region containing the PBH and the fiat Friedmann region is 
made inside the cosmological particle horizon of the flat Friedmann solution, then we describe the PBH as "causally 
produced" . 

If the perturbation is produced by some hydrodynamical process, it is natural to assume that the perturbation is 
confined within the cosmological sonic horizon, which is the distance traversed by a sound-wave since the big bang. 
In a fiat Friedmann universe containing a perfect fluid with the equation of state p = fcp (0 < fc < 1), this is just ^/k 
times the particle horizon size. For a stiff equation of state A: = 1 or a massless scalar field, the cosmological sonic 
horizon coincides with the particle horizon. 

Definition 3 // the matching between the perturbed region containing the PBH and the flat Friedmann region is made 
inside the cosmological sonic horizon of the flat Friedmann solution, then we describe the PBH as "hydrodynamically 
produced" . 
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In fact, this is a very natural value for the matching radius, since the existence of a sound-wave allows for a discontinuity 
in the pressure and density gradients 0. Note that the cosmological sonic horizon is somewhat smaller than the 
Jeans length, which is the scale on which perturbations oscillate as sound-waves rather than grow. As discussed in 
the Introduction, an overdense region needs to be bigger than the Jeans length to collapse against the pressure 
and this might seem to be incompatible with the last definition. However, in this argument, the relevant Jeans length 
is measured at the time when the overdense region has its maximum expansion and this is somewhat before the black 
hole itself forms. 

We now calculate these scales explicitly for a flat Friedmann universe with 

ds^ = -dt^ + a'^{t)[dx^ + X^{d0^ + si-a^ 9d(l)^)]. (4.1) 

In this case, the null coordinates are related to the conformal time and comoving radial coordinates i] and x by 

u = Tj-x, (4.2) 
V = V + X, (4-3) 

where adrj = dt. It is most convenient to derive expressions in terms of u and v since this is more suitable for the 
application of the results obtained in the previous section. In such coordinates the flat Friedmann spacetime for a 
perfect fluid with p = kp is given by 

/,,,,, \2/(l+3fc) 

a = Crj^/i^+m = c(^] , (4.4) 



/ , X 2/(l+3fe) 
„ I U + V\ " ' V - U 

where C(> 0) is a constant of integration. In the following we assume < < 1. 
The cosmological particle horizon is given by 

ucMv) = 0, (4.6) 
(i.e. r] = x)- Using Eas. H2.10|) . H4.4|l and H4.5|l . we obtain 

2m / 2 , „ „ , , 

(4.7) 



so the cosmological apparent horizon has the retarded time coordinate 



ucAHiv) = - I ) (4-8) 



1 



Therefore, the cosmological apparent horizon is spacelike for 1/3 < /c < 1, null for k — 1/3 and timelike for < < 1/3, 
as seen in Fig. El For < fc < 1, the cosmological sonic horizon has the retarded coordinate 



v-x 



1-Vk 



ucsh(w) = — ; — V = jz^v, (4.9) 

\V + Xj 1 + Vk 

since the sound-speed is yjdp/dp ~ \fk. As seen in Fig. |21 the cosmological sonic horizon is timelike and inside the 
cosmological particle horizon for < fc < 1, while it is null and coincides with the cosmological particle horizon for 
A; = 1. It is also null for a flat Friedmann spacetime with a massless scalar field, so ucsh = in these cases. A similar 
figure has been presented in [20| . 

The area radii for these three cosmological horizons are written in terms of the advanced time coordinate v as 

rcPH(^) = CQ)™, (4.10) 



rCAH(w) = C 



l4-3fc 



2 

C\/fc I - 



2v 



3(1 + fc) 

3(l + fc) 



(4.11) 



rcsH(«) = G^k^j—j=\ , (4.12) 



6 



where we use the abbreviation rcAniv) = r(ucAH(w), v). It is also convenient to express these scales in terms of the 
Hubble parameter H. Since this is defined as 

H.Y-^.-l-±. (413, 

we have 

aH=- -^7 (4-14) 

The area radii of the cosmological particle horizon, apparent horizon and sonic horizon are then given by 

rcAH(«) = H-\ (4.16) 

rcsKiv) = ^^H-' (4.17) 

in a flat Friedmann universe. Note that H is here the Hubble parameter at the cosmological apparent horizon of the 
exact Friedmann background and this may not be the same as the Hubble parameter within the perturbed region 
itself. 

Another important quantity is the Jeans scale, which is given by |2l| 



Aj = W^, 4.18 

\l P + P 

although this expression is strictly only valid for Aj less than the cosmological apparent horizon size. Below this 
scale density fluctuations oscillate as sound-waves rather than growing. However, it should be stressed that this is 
the Jeans length in the background universe, which is different from the Jeans length in an overdense region. It is 
usually assumed that a region can collapse to a black hole only if is larger than Jeans length at maximum expansion 
but, as discussed later, this is somewhat less than the value given by above. Since Aj is a wavelength, the radius 
corresponding to Eq. (|4.18|) is 



1 / 2k 



where we have used the Friedmann relation 

= fp. (4.20) 

However, we note that Coles |2^ gives a different equation for the Jeans length in a fluid with equation of state 
parameter k and this leads to the expression 

which we henceforth adopt. This is always of order the sonic horizon but it is not exactly the same. More precisely, 
the ratio of the physical scales given by Eqs. (|4.17|) and (|4.21|) is 

^ ^ 5 + 9fc 
rj 27r(l + 3fc) ^ ' 

As k goes from to 1, this decreases monotonically from 0.796 to 0.477. We should note that the nature of the Jeans 
instability will be considerably modified for k > 0.1 because of the cosmological expansion 21]. These scales are 
depicted in Fig.|3| It is interesting to note that rcpH — fcAH at fc = 1/3 and rcsH has a maximum there. 
At the cosmological particle horizon, equation (|4.7(l with u = Q implies 



2m 



r \l + 3k 



(4.23) 
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Therefore 2m/r < 1 for 1/3 < < 1, 2m/r = 1 for k = 1/3 and 2m/r > 1 for < A: < 1/3. This means that, for 
1/3 < fc < 1, the cosmological particle horizon is inside the cosmofogical apparent horizon, i.e. ucph > ^cah- In this 
case, r > 0, < and r,^ > at the cosmofogical particle horizon. Since a causally-produced PBH is associated 
with a perturbed region inside the cosmological particle horizon. Theorem |31 immediately implies the following: 

Theorem 4 For a PBH causally produced from a perfect fluid with p > and p = kp (1/3 < k < 1), the ratio of the 
area radii of the black hole and cosmological apparent horizons on the same ingoing null surface is less than 



( ' V 


'3(1 + A:)" 


Vi + 3fcy 
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Proof. For 1/3 < /c < 1, we have r > 0, r^^, > and r.„ < at (ucph('I'), v). Theorem then gives 

. N ^ o / N 2mcpH(wo) rcpniv) 



2 



1 + 3/c 



l + 3k 



3(1 + fc) 



where we have used Eqs. (|4.10ll and (|4.11() . □ 



This upper limit is plotted in Fig. O and is always less than the cosmological particle horizon for k > 1/3. Note that 
the ratio in Theorem Speaks with a value of 1 for a radiation equation of state and this is the most natural one for 
the early Universe. 

Theorem 5 For a PBH causally produced from a massless scalar field, the ratio of the area radii of the black hole 
and cosmological apparent horizons on the same ingoing null surface is less than 3-\/6/32. 

Proof. This follows trivially by putting fc = 1 in Theorem 0] □ 



For 1/3 < fc < 1 or a massless scalar field, once we have fixed the matching coordinate Wmatch, the radius of the 
black hole apparent horizon is maximized if there is a vacuum outside this horizon. If we also change Umatch, the 
upper limit on the black hole size is attained by the vacuum solution in the limit Umatch — > wcph- For < fc < 1/3, 
the radius of a causally produced PBH could be as large as the cosmological particle horizon size. 

At the cosmological sonic horizon, we have 



2m 
r 



1 + 3fc 



(4.24) 



and this is necessarily less than 1 for < fc < 1. In this case, r „ < and r ^, > at the cosmological sonic 
horizon. For fc = 1, the cosmological sonic horizon coincides with the cosmological particle horizon. Therefore, for a 
hydrodynamically-produced PBH, we can prove the following. 

Theorem 6 For a PBH produced hydrodynamically from a perfect fluid with p > and p — kp (0 < k < 1), the ratio 
of the area radii of the black hole and cosmological apparent horizons on the same ingoing null surface is less than: 



( ' V 


r 3(1 


4-fc) " 


Vi + 3fcy 


.2(14 


-\/fc)J 



-| 3(l+fc)/(l+3fe) 

Proof. For < fc < 1, we have r > 0, > and < at (wcsh(w), v). Then we have 

2mcsH(w) rcsniv) 



(4.25) 



rBuiv) < 2mcsH(w) 



fcSH(w) 7'cah(w) 



fCAH(w) 



l + 3fc 

where we have used Eqs. (|4.11|l and (|4.12l) . □ 



' 2Vk ' 
1 + 3fc 



3(1 + fc) 



2(1 + Vfc) 



3(l+fe)/(l+3fc) 



fCAH(w), 



This upper limit is also plotted in Fig. O It is noted that the upper limit (|4.25|) takes a maximum value of 0.309 at 
fc ~ 0.343. 
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V. NON-LOCALISED PERTURBATIONS AND SEPARATE UNIVERSE CONDITION 



We have considered PBHs formed from localised perturbations. However, since the Friedmann universe will be 
perturbed everywhere in some situations (e.g. if the perturbations derive from inflation), the above upper limit 
analysis may not always apply. In this section we derive the condition that a PBH is part of our universe rather than 
being a separate closed universe. This argument does not depend on whether the perturbation is localised or not. 
For clarity we include factors c and G explicitly. To make the discussion more precise, we assume that the overdense 
region which evolves to the PBH is part of a homogeneous closed Friedmann model with 

ds^ = -c^df + a^{t)[dx^ + sin^ x{d9^ + sin^ Ocf)^)]. (5.1) 

The cosmic scale factor within this region then evolves according to 

'd\^ SnGp (? 



(5.2) 



where a dot denotes the derivative with respect to t. This implies that the scale factor at maximum expansion is 

(5.3) 



For an equation of state parameter A:, the background Friedmann density is 

If we denote the overdensity at maximum expansion by A = (pmax/pb), equation H5.3|l then gives 

a„,ax = ^(1 + fc)A-i/2ci = A-i/2ci7-i. (5.5) 

Since the maximum separation between two points on a spatial hypersphere is tt times the radius of the hypersphere, 
the largest possible radius for a region at maximum expansion is rmax — T^o-n 
The overdensity parameter A is expressed in terms of k as 



^max ■ 



A = ^ 

4 



3(1 + fc) 



1 + 3/c 



2 r 



{ 3(i+fc) y 

\^2(l+3fc) J 
^ \l+3k J 



(5.6) 



where the derivation is given in Appendix. From Eq. H5.5|l our final expression for the separate universe scale in units 
of iJ-i is 

= = 2V^ (l + 3fc) li+3. j 

^ 3(1 + k)rf ^ ^ ' 

1,2(1+3/=) ) 

In a dust universe (fc = 0), this gives 4/3. In a radiation universe (fc = 1/3), it gives 7r/2. In an (unphysical) universe 
with fc = oo, it gives 2. This expression is shown in Figure El Note that rmin exceeds the cosmological particle 
horizon for k > 0.126. It should be stressed that, although we have presented our analysis in terms of the evolution 
of the overdense region, such a region does not evolve into separate universe. If it is a separate universe at maximum 
expansion, it is always separate. 

We have noted that the Jeans length within an overdense region at maximum expansion is less than the background 
Jeans length. Since the Jeans length depends upon the inverse square root of the density, one might expect the 
reduction to be of order A^/^, although this is not an exact result since the expression given by Eq. H4.18|l only applies 
in the linear regime. The important point is that the usual lower limit on the size of a black hole at formation is likely 
to be below the upper limits shown in Fig. |3 Furthermore the existence of critical phenomena may circumvent the 
Jeans limit anyway, since this permits black holes of arbitrarily small mass. 
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VI. PBH ACCRETION 

If we neglect the background expansion and assume a quasi-stationary flow, a simple Newtonian treatment for 
a general fluid suggests that the accretion rate for a black hole of mass M should be given approximately by 

— =4TrapRlvs, (6.1) 

where Ra = GM /v^ is the accretion radius, Vs is the sound-speed in the accreted fluid, a is a factor of order unity 
which depends upon other physical features of the flow, and factors G and c are now written explicitly. For a flat 
Friedmann universe containing fluid with equation of state p — kp, we have 



Vkc, (6.2) 
GAj 
kc^ 



Ra = (6.3) 



Therefore the accretion rate is 



This can be integrated to give 



67rG(l + kff^ ' 
dM _ 2aGM^ 

~ 3fc3/2(l + fc)2c3t2 



(6.4) 



(6.5) 



At 

M = ^, (6.6) 

A . ^^^^!(i±^, (6.7) 
2aG 

where M{ is the black hole mass at the time ti of formation. If we assume that Mf = ^At{ with ^ < 1, then Eq. H6.6|l 
implies 

M^Mf(l-0~^ as t^oo. (6.8) 

If ^ ^ 1, the black hole could not grow very much. However, if ^ is close to 1, then the black hole could grow 
significantly. In other words, if the black hole mass is initially very close to At{, which is approximately the same as 
the mass contained within the cosmological sonic horizon, the black hole can grow considerably. In particular, in the 
limit ^ = f , Eq. H6.6|l implies M c>c so the black hole grows at the same rate as the universe. The physical scale 
corresponding to the mass for which this analysis predicts self-similar growth is given by 

rfff = 2a"^(H-fc)fc3/2. (6.9) 

This scale is also indicated in Fig.|21(for a = 1) and one sees that in this case it violates the separate universe condition 
for sufficiently large values of k. This is another reason for regarding the Newtonian prediction as suspect. 

Since the above calculation neglects the effects of the cosmological expansion, one needs a relativistic calculation 
to check this. The Newtonian result suggests that one should look for a spherically symmetric self-similar solution, 
in which every dimensionless variable is a function of z = r/t, so that it is unchanged by the transformation t — > 
at, r ar for any constant a. This problem has an interesting but rather convolved history. By looking for a black 
hole solution attached to an exact Friedmann solution via a sonic point, Carr & Hawking first showed that there is 
no such solution for a radiation fluid ^] and the argument was later extended jlfl] to a general p = kp fluid with 
< k < 1. Lin et al. M subsequently claimed that there is such a solution in the special case fc = 1. However, 
Bicknell & Hcnriksen uz| then showed that this solution is unphysical, in that the density gradient diverges at the 
event horizon. The solution can be completed only by attaching a Vaidya ingoing radiation solution interior to some 
surface (i.e. the scalar field has to turn into a null fluid). This is rather contrived, which suggests that the black 
hole must soon become much smaller than the particle horizon, after which Eq. H6.1|l implies there will be very little 
further accretion. 

The present results are qualitatively consistent with numerical calculations [T^ Il4l [Tsl Il6j. In particular, simu- 
lations 01 suggest that, the harder the equation of state, the more difficult it will be for a PBH to form and the 
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smaller the mass accretion rate will be. This may be understood in the context of the present analysis. Even if the 
initial perturbation is larger than the cosmological apparent horizon, as may arise in the inflationary scenario, it soon 
becomes comparable with the cosmological particle horizon, since the propagation speed of the perturbation is the 
sound speed, and this is slower than the light speed if A; < 1. This means that the constraint on causally produced 
PBHs obtained also applies in this case. We conclude that the maximum mass and radius of these PBHs will decrease 
as the equation of state becomes harder if fc > 1/3. 

VII. CONCLUSION 

We have proved that a causally or hydrodynamically produced PBH cannot be as large as the cosmological apparent 
horizon for either a perfect fluid with a hard equation of state or a massless scalar field. It is clear that the result can 
be extended to a wider class of matter fields or to more general cosmological solutions. We have also derived an exact 
expression for the maximum size that an overdense region can have without being a separate universe. This is always 
larger than the cosmological apparent horizon and larger than the particle horizon size for a sufficiently hard equation 
of state. Finally, we have considered accretion by PBHs but numerical simulations will be necessary to determine this 
precisely. For the perfect fiuid case, such simulations indicate a small accretion rate. For a massless scalar field, work 
is in progress by the present authors 0. 

APPENDIX: CALCULATION OF THE OVERDENSITY A 

We here calculate the parameter A appearing in Eq. H5.5|) . In the dust case this is well known to be (37r/4)^ but 
we need to generalize this result. For general fc, we have p oc a"^^^"'"'^-* and so Eq. I|5.2|l can be written as 

a2 ^ ^a-(i+3'=) _ B (A.l) 

where 

A = d^^n,al+^\ B = al{n,-l)^^. (A.2) 

Here the subscript "□" indicates some initial epoch at which the overdensity is small, so that the density parameter 
is close to unity. 

By defining a new scale factor b and a new time coordinate r such that 

6 = 0^+3^, dr ^ {l + ?,k)b^dt, (A.3) 
we can transform Eq. ljA.l|) into the dust form: 



dhV A 



dr 



--c\ (A.4) 



This has the well-known parametric solution: 



b = ^&max(l - COS 77), T = -r,„ax(?7 " Sm T]) (A. 5) 

Z TT 

where rj = n corresponds to the epoch of maximum expansion. Equations (jA.3|l and (jA.4p then imply 

t = / b i+'if^dr 



1 + 3fc 
1 

1 + 3fc 







(Il^^j J (l-cos7y)TWd,7. (A.6) 



Equations l|A.l|l and (|A.2|I give 
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so we have 



-'max / ^max 



Po 



-3(l+fc) 



3(l + fc) 
f^n - 1\ ' + 



Since Eq. (|5.4() implies that the background density is given by 



Ph 



Po \ ^max , 

where to and imax are given by Eq. HA.6|) , the overdensity can be expressed as 



A = 



/J'(l - COS 77) 1+3'' d?] 

( 1 — COS ?/) TTsfc c/r/ 



Equation (|A.5p impHes that the integral limit rjo, corresponding to the time to, is given by 



77o = cos 

The top integral in Eq. IjA.lOp can be expressed as 



/■TT 

/ (1 - cos77)TT3fcd?7 = 2TT3fc ^ 
^0 



»o- 1 



p /i(l±fc) \ 
^ \^2(i+3fc) y 

p 2+3fc\ 



(A., 



(A.9) 



(A.IO) 



{A.11) 



(A.12) 



where T{z) denotes the gamma function. By using the fact that rfa is small, the lower integral in Eq. (|A.10p can be 
simplified to 



3(l + fc) 



3(l + /i:)y V 



(A.13) 



where we have used Eq. IjA.lip . One observes that the terms involving Vlo in the expression for A necessarily cancel 
to give 



A 



3(1 + fc) 



i2 r 



1 + 3fc 



( 3(i+fc) y 

\^2(l+3fc) ) 

Y ( 2+3fc\^ 
^ \^ l+3fc ) 



(A. 14) 
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FIG. 1: Schematic figure for proofs of Theorems and Lemma|21 Here Uc and mbh denote the centre and the black hole 
apparent horizon, respectively. 
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FIG. 2: The coiiformal diagrams for flat Fricdmann spacctimcs with a perfect fluid with tfie equation of state p = kp for (a) 
< fc < 1/3, (b) fc = 1/3, (c) 1/3 < fc < 1 and (d) fc = 1 or a massless scalar field. CAH, CPH, CSH denote the cosmological 
apparent horizon, particle horizon and sonic horizon, respectively. 
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FIG. 3: The characteristic physical scales in units oi for different values of k, where H is the Hubble parameter at the 
cosmological apparent horizon of the exact Friedmann background. CPH, CAH and CSH denote the cosmological particle 
horizon, apparent horizon and sonic horizon, respectively. "Jeans" denotes the Jeans length implied by Coles |23I . These four 
scales are plotted as thin lines. Causal and Hydro denote the upper limits for causally and hydrodynamically produced PBHs, 
respectively. SU denotes the upper limit which comes from the separate universe condition. These three scales are plotted as 
thick lines. SS denotes the black hole size which the Newtonian argument of self-similar growth suggests. 



